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QQ Abstract 

I In this paper in the space L'^\o, 1) the problem of construction of optimal quadrature formulas is 

CNJ , considered. Here the quadrature sum consists on values of integrand at nodes and values of first derivative 

of integrand at the end points of integration interval. The optimal coefficients are found and norm of the 
, error functional is calculated for arbitrary fixed A'^ and for any m > 2. It is shown that when m = 2 and 

m — 3 the Euler-Maclaurin quadrature formula is optimal. 
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It is known, that numerical integration formulae, or quadrature formulae, are methods for the 

approximate evaluation of definite integrals. They are needed for the computation of those 

^ . integrals for which either the antiderivative of the integrand cannot be expressed in terms of 

(N 

elementary functions or for which the integrand is available only at discrete points, for example 
OO ! 

O ■ from experimental data. In addition and even more important, quadrature formulae provide a 

- ^ . basic and important tool for the numerical solution of differential and integral equations. 

H ' There are various methods in the theory of quadrature, which allow us approximately calcu- 

late integrals with the help of finite number values of integrand. Present work also is devoted 
to one of such methods, i.e. to construction of optimal quadrature formulas for approximate 
evaluation of definite integrals in the space L'^\o, 1) equipped with the norm 
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Consider following quadrature formula 



^ N 

I ^{x)dx=Y,C[mP] + A^'[Q]+B^'[N] (1.1) 

n /3=0 



1 



with the error functional 



N 



e(x) = £[0,1] (x) - - hp) + A6'{x) + B6'{x - 1) 

(3=0 



[1.2) 



in the space L2"^(0, 1). Here C[/3], P = 0, N, A and B are the coefficients of the formula (1.1), 
[P] — hp, h — j^, N — 1,2,... , £[0,1] (x) is the indicator of interval [0,1], S{x) is the Dirac 
delta-function. 
The difference 

\ 

{i{x), ^(x)) = / ^{x)dx- C[PHP] - V[0] - B^'[N]. 
/^=o 
is called the error of the quadrature formula (1.1) 

Error of the formula (1.1) is estimated with the help of norm of the error functional (1.2) in 
the conjugate space L^^*{0, 1), i.e. by 



e{x)\L^r^* 



sup \{e{x),ip{x))\. 



\mx)\L. 



(m) 



Furthermore, norm of the error functional £{x) depends on the coefficients C[/3], A and B. 
Choice of the coefficients when nodes are fixed is linear problem. Therefore we minimize norm 
of the functional i{x) by coefficients, i.e. we find 



e{x)\L 



(m)* 
2 



= inf 

Cll3],A,B 



e{x)\L 



(rn)* 
2 
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If 



£{x)\L2' is found then the functional £{x) is said to be correspond to the optimal quadra- 
ture formula (1.1) in l'^^ and corresponding coefficients are called optimal. Thus we get fol- 
lowing problems. 

Problem 1. Find norm of the error functional £{x) of quadrature formula of the form (1-1) 

in the space L^^'^*{'^), 1). 

Problem 2. Find coefficients C[/9], A and B which satisfy the equality (1-3). 
Problem 2 for quadrature formulas of the form 



N 



N 



'x)dx = y^^Pk(p{k) 



k=0 



on Ly' first considered by A.Sard [1]. By A.Sard and S.D.Meyers [2] the solution of this 
problem was obtained for the following cases: m = 1 for arbitrary fixed N; m — 2 ior N < 20; 
m = 3 for iV < 12; m = 4 for < 9. 

By I.J.Shoenberg and S.D.Silliman [3] the Sard's problem for N oo, i.e. for formula of 
the form 



is considered. In [3] an algorithm for finding of the coefficients B^^' is given with the help of 
spline of degree 2m — 1. In the cases m = 2, 3, 7 the coefficients S^™^ are calculated using a 
Computer. Calculation of these coefficients up to m = 30 were done by F.Ya.Zagirova [4]. 
In [5] in the space l'^^ considered quadrature formula of the form 



where < r^j < 1 , u{x) is weight function, C[P] are the coefficients. In [5] the algorithm for 
finding optimal coefficients C[/3] of quadrature formulas of the form (1.4) is given and for the op- 
timal coefficients the system of 2m — 2 hnear equations is obtained. These results of S.L.Sobolev 
generalized above mentioned results of A.Sard, S.D.Meyers, Schoenberg and Silliman. Further 
realization of Sobolev's algorithm studied by Z.Jamalov, F.Ya.Zagirova, Kh.M.Shadimetov. In 
[6], [7] the problem of construction of optimal quadrature formulas (1.4) was completely solved 
for arbitrarily fixed N and for any m in the space l'^\o, 1). 



Main goal of the present work is to solve problems 1 and 2 for quadrature formulas of the 
form (1.1). 



2 Definitions and known formulas 

In this section we give some definitions and formulas which are necessary in the proofs of main 
results. 




k=0 








(1.4) 
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Euler polynomials E}.{x) , A; = 1, 2, ... is defined by following formula [8] 

Eo{x) = 1. 

For Euler polynomials following identity hold 

Ek{x) = x'^Ek (^i^ , (2.2) 

and also following theorem is take placed 

Theorem 2.1 [9]. Polynomial Pk{x) which determined by formula 

= (2.3) 

is the Euler polynomial (2.1) of degree k, i.e. Pk{x) — Ek{x). 
Following formula is valid [10]: 

E - E A'O' - E (2.4) 

7=0 ^ 1=0 V ^ / ^ i=o V ^ / 

where A*7'^ is finite difference of order i of 7*^, A^O*^ = A*7'^|^=o- 
At last we give following well known formulas from [11] 

^-1 fc+i , , „ 

E^'-Ej^ff^'^'. (^-s) 

7=0 j=l J V ' J J 

where B^+i^j are Bernoulli numbers, 



p=o\P J 

3 The extremal function and the representation of the error func- 
tional norm 

To solve problem 1, i.e. for finding norm of the error functional (1.2) in the space L^\o, 1) 
concept of the extremal function is used [12]. The function ■ip£{x) is said to be extremal function 



of the error functional (1.2) (see [12]) if following equality holds 



(m) 



In the space L^\o,l) the extremal function V'^(x) of the error functional £{x) is found by 
S.L.Sobolev. This extremal function have the form 



M^) = i-^rK^) * G{X) + Pm-l{x), 



(3.1) 



where G{x) 



signx 



2(2m-i)! ' ^"1-1(3:^) is a polynomial of degree m — 1. Since the functional l{x) 
belongs to the space L^'^-**(0, 1) therefore following holds 



(£(a;),x") = 0, a = 0, 1, m - 1. 



(3.2) 



Norm of the error functional of quadrature formula (1.1) depends on coefficients of this formula. 

Indeed, since the space 4'"^(0, 1) is Hilbcrt space, then by using (3.1), taking into account of 
Riesz theorem about common form of a linear continuous functional on Hilbert space, we get 



= (£,^,) = (-l) 



m+l 



A-B 



2 A 



(2m -3)! I J 2(2m-2)! 





X 



2m-2 



dx- B 



1 

/ 



{x-l) 



2m- 2 



2(2m-2)! 



■dx + 



N 



+2Y,m[A 



/3=0 



2(2m-2)! 2(2m-2)! 



N ^ 



\x - h(3\ 



2m- 1 



2(2m- 1)! 



-dx- 



/3=0 7=0 ^ ^00 



X — yY'^ ^sign{x — y) 
2(2m- 1)! 



dxdy 



Thus, the problem 1 is solved for quadrature formulas of the form (1.1) in the space L2™''(0, 1). 



4 The system of Wiener- Hopf type 

Now we investigate problem 2. For finding of minimum of the ||£||^ under the conditions (3.2) 
Lagrange method of undetermined multipliers is used. For this we consider following function 

m— 1 

* = ||£||^ + 2-(-l)-+i J]A„(£,a;"). 
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Equating to zero partial derivatives by coefficients C[(3], A and B, togetfier witfi conditions 
(3.2) we get following system of linear equations 



2(2rn-l)! 2(2m-2)! 



^*^^^^2(2m-2)! ^2(2m-3)! ~ " 2(2m-l)!' ^"^'^^ 



E - 2(2m-3)! + £ ^ 2(2m-l)! ' ^'-'^ 

7=0 

^ 1 

J]C[7]/i7 + ^ + 5=^, (4.5) 



7=0 \ / \ / a=l 



7=0 

N 



VC[7](/i7)- + «5=^, Q; = 2,m-1. (4.6) 
The system (4.1)-(4.6) is called by system of Wiener- Hof type for the optimal coefficients [12]. 



In the system (4.1)-(4.6) coefficients C[P], P — 0, N, A and B, and also Aq, a — 0,m — 1 are 
unknowns. The system (4.1)-(4.6) has unique solution. The proof of existence and uniqueness 
of the solution of this system is as the proof of existence and uniqueness of the solution of 
Wiener-Hopf system of the optimal coefficients in the space L'^\o, 1) for quadrature formulas 
of the form (1.4) (see [13]). 



5 The optimal coefficients and norm of the error functional 

In present section we study solution of the system (4.1)-(4.6). In the solution of this system 
we use the approach which used in solution of the linear system for optimal coefficients of 
quadrature formulas of the form (1.4) in [6]. 
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5.1 The coefficients of optimal quadrature formulas 

It is easy to prove following lemma for the coefficients C[/3] of quadrature formulas of the form 
(1.1). 

Lemma 5.1. The optimal coefficients C[P], 1 < (3 < N — 1, of quadrature formulas of the 
form (1.1) have following form 

(m-l \ 
i+^(4gf+Pfegr^)j, i</?<iv-i, (5.1) 

where dk,Pk are unknowns, are roots of the Euler polynomial E2m-2{<l), llkl < 1- 

Lemma is proved as lemma 3 of the work [9] and in the proof the discrete analogue Dyn[(3] 

of the polyharmonic operator is used. The discrete analogue Dm[P] of the polyharmonic 

operator is constructed in [14]. 

We need following lemmas in proof of main results. 
Lemma 5.2. Following identity is take placed 

± ^^7^'':''-''"'a-0° ^ (-1)-. ± ''^'+"~"'-^)'"a-0°. (6.2) 

^ (q- 1)^+1 ^ ^ ^ (1 - g)'+i ^ ' 

i=o > j=0 ^ ^' 

here a and N are natural numbers, A*0" is finite difference of order i of 7" at the point 0. 

Proof. For convenience left and right sides of (5.2)we denote by Li L2 respectively, i.e. 

_^ dq + pq^^^{-iy^\ ,.a^i^ dq^+pq^^\-iy^\ 

^i-Z. ^-Tp ^° ^2-(-i) 2. AO. 

i=o ' i=o ^ 

First consider Li. By using the equality (2.3) and identity (2.2) for Li consequently we get 

^ dq + pq^+^{-iy+\ , dq ^^+"(-1)"+^ /l\ 

= L (g-l)^^^ ^ ' - (^31)^^"-^'^) + (.-1)-^ """-^ UJ ^ 



dq pg^+"(-l)"+^£;,_i(g) rfg + pg^+^(-l)°+^ ^ 

+ — T^T^Ti 3^=^ = 7- TTT^T^ Ec,-i{q). (5.3) 



(g_l)a+i "-^vi/ (g-l)«+i ga-i (g-l)"+i 
Similarly for L2 by using (2.3) and (2.2) we have 
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-Ea-l{q) + TZ TT— T^a-l(g) = 7- TT—. ^a-l(g) 



From (5.3) and (5.4) clear, that Li = (— 1)"+^L2. Lemma 5.2 is proved. 
We denote 



Lemma 5.3. Following identities are valid 

m—l I — 1 2m— 3- 



and 



^7r£^zr(2;;rT^^7^+ .£;^7r ^ i\ (2m-i-j-z)! 

^ (?■ - 1)! ^ d 



p^(i-l)! ^ p! (2m- 1 - J -p)! 

The proof of lemma is obtained by expansion in powers of h of left sides of given identities. 

For the coefficients of optimal quadrature formulas of the form (1.1) following theorem holds. 
Theorem 5.1. Among quadrature formulas of the form (1-1) with the error functional (1-2) 
there exists unique optimal formula which coefficients are determined by following formulas 



pi°i-''(^+|:' '"f_",f'' ). (8-6) 

m-l \ 

C[(3] ^h[l + J2 {dkqi + Pkqu'^) , /? = T;A^, (5.6) 



fe=i 

m-l , 



k=l ^'^ 



^-^[l2 ^ (1-,.)^ j' ^'-'^ 



n-h^ 1 I V- ^fcgfc +Pkqk \ 



where and Pk satisfy following system 2m — 2 linear equations: 

+ , = (5.10) 



m-12m-2 , , -,Af+«/'_1 V+1 



E E '''""•r''!\vi7"' ^'"'""^ = ». (5-11) 

fc=l 1=0 ^^'^ ^ 

E E (1 - '^") ^"HT- fv^i"'''' ^'"^' - °' (5.12) 
fe=i i=o 

m-l 2m-2 , ^ , jv+j , „ „ 
fc=l 1=0 ^ 

i/ere are Bernoulli numbers, A^'-f^ is difference of order i of A*0-^ = A*7-' |-y=o , qk are 
roots of Euler polynomial of degree 2m — 2 ,\qk\ < 1. 
Proof. First we give plan of proof. 



Prom (5.1) clear that instead of unknowns C[P], /? = 1, A?" — 1 it is sufficient to find unknowns 



dk, Pk, k — l,m — 1. The coefficients C[0], C[N], A, B and A^, o; = 0, m — 1 are expressed by 



dk and pk-, k ~ l^m — 1. So if we find dk and pk-, then the system (4.1)-(4.6) is solved completely. 
Substituting the equality (5.1) to equation (4.1) we get polynomial of degree 2m of hj3 on both 



sides of (4.1). Equating coefficients of same degrees of hjd we find Aq,, a = 0, m — 1, C[0], A and 
system (5.10) for dk, Pk- Taking account of (5.1), (5.5), (5.9), from conditions (4.4) and (4.5) 
we get (5.7), (5.9), i.e. we obtain C[N\ and B. Further, by using (5.1), (5.9) and expression 
for Ai, from (4.2) we get the equation (5.11). System of equations (5.12) for unknowns dk ,Pk , 
we obtain from equation (4.6), using (5.1), (5.5)-(5.9). Finally, taking into account (5.1), (5.8) 



and Xa, a = l, m — 1 from equation (4.3) we have (5.13). 

Further we give detailed explanation of proof of the theorem. 
First we consider first sum of equation (4.1). For this sum we have 

^-E^W 2(2m-l)! ^ 
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^^(2;^ + (2m -1)! -^^W 2(2m-l)! 



Lat two sums of the expression S we denote 



/9 



^^^S'^^"^ (2—1)! ' = 2(2^-1). 



7=i 7: 

and we calculate them separately. 

By using lemma 5.1 and formulas (2.4), (2.5) for Si we have 

k=l 



p I m— i 



7=0 



(2m- 1)! 



^2m 



(2m- 1)! 

//-'"' 



/3-1 m-1 / /3-1 

E 7^"^-^ + E ( E ^rf "^-^ + vk€-' E ^.'7'"^"' 

7=0 7=0 



.7=0 

2m 



ifc=l 



(2m- 1)! 



^ j!-(2m-j)! ^""^ 



m— 1 



2m- 1 



, ^ J 'k ATI"-' 



1-/3 2m-l , j^2m-l 
-1 (^^-1)' 



N-8 

PkQk 



1 



2m- 1 



2m- 1 



^2m— 1 



9fc - 1 



2m- 1 



Taking into account that is a root of Eulaer polynomial E2m-2{'l) an using formulas (2.3), 
(2.6) the expression for we reduce to following form 

\2m rhn\2m-l 2m-2 



Si 



m 



(2m)! (2m-l)! ' ^j!(2m-j)! 



2m- 1 



o2m-l-j 3 rl n J- ^ rr^+'i-f 1 
il('9m _ 1 _ 7)! ' 



j\(2m-l-j)\ 



(5.14) 



(Qk - iy+' 

Now consider 5*2. By using conditions of orthogonality (4.4)-(4.6) the expression S2 we rewrite 
by powers of h(3 

N 



7=0 



2m- 1 ^ '^-^ ^f^pym-l-j 



E 



{hf^ - ^7) 

2(2m-l)! 2^j!(2m-l-i)! + l 



-jB]- 
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ahw. [-2-^-^)^ Uhry. ^ 2 E ,,('J_i_,)| EgM(-M^ (5.15, 



The right side of the equation (4.1) have following form 



J 2(2ro-l ! (2ro ! ^ 2(2m - 1 - j)!y + 1 ! ^ ' 

Substituting (5.16) and 5" into equation (4.1) and using (5.14), (5.15) we have 

(2m)! ^ ^(2m-l)! (2m-l)! ' j\{2m-j)\ 
^ j!(2m-l-j)! trt^ 

_ 1 (/./j)^-^-^(-i)^- / 1 _ \ m^-' (i-A-R 

2 2^ j\{2m-l-j)\ + 1 ^ ;^2(2m-2)!V2 

J— ^ 

2(2m-l)! 2 j!(2m-l-i)! '^^ 2(2m-2)! 

(fe/3)^-^-^(-l)^- , (-/i/?)^-^-^- 

^ 2 . j! . (2m - 2 - i)! + ^ '''^^'^^^ (2m)! + A. 2 • (2m - 1 - j)! • (j + 1)! 

Hence equating coefficients of same powers of /i/? gives 

, _ 1 ( {-If^-^ B2^-jh^"'-^ 



^ {2m - I - j)\ ■ j\ \2{2m - j) 2m - j 



2. (g, - 1)^+1 ^0 + 

fc=l i=0 ^^'^ ^ 

+ 2 E CbK-hjY'^-'-^ _ (2m - 1 - j) • 5 • j , J = 1, 2, m - 1, (5.17) 

7=0 ^ 

11^ 1 

^°=2T(SS^ + 2T(2;;73T^gC'b](-M--'-B^:p^ (5.18) 

trs (*-i)'+' j + 1 
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Here equation (5.19) is the equation (5.10) for unknowns dk and pk- 

Substituting expressions (5.20) and (5.21) into (4.4) and (4.5), also taking into account (5.1), 
we find C[N\ and B, wfiicli fiave following form 

(l^E^^^ff^), (5.22, 

Now substituting the expression of Ai from (5.17) when j — 1 into (4.2), we get one more 
equation with respect to unknowns dk and pk- 

k=l i=0 ^^'^ > 

i.e. we obtain the equation (5.11). 

Next, to obtain (5.12) we use equation (4.6). Consider equation (4.6). Since a = 2,m — 1, 
then 

We denote L = X) C[j]{hlT- Using (5.1), (2.4), (2.5) for L we get 

7=1 

(7V-1 m-1 / AT-l 7V-1 

E + E E + ^'^^^'^ E ^^7" 

7=1 fc=l \ 7=1 7=1 



^a+1 QfcQfc +pfcgfc(-ij ^^^^^^ 



fe=l i=0 
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Hence taking into account(2.6) and grouping in powers of h, we have 



Substitution of obtained expression of L to (5.25) gives 



a + l^^(7-l)! (a + l-?)! I 7 ^ 



' m—l a 



Hence keeping in mind (5.22), (5.23) we get 



2-.(j-l)\(a + l-j)\[j (1-9.)^+^ 



+^""' (E E (1 - ) ^ '|^",'l'i^.+/''' ^^0") = 0- (5-26) 

\ /c^l i — / 

Clearly that the left side of (5.26) is polynomial of degree a + 1 with respect to h. Prom (5.26) 
we obtain that each coefficient of this polynomial is equal to zero, i.e. 



k=l i=0 

and 



EE (1 - gf) ^ 'd'i+i''''' ^'0" = 0- (5-28) 

fe=l i=o ^> 

Since a = 2, m — 1 then from (5.27) and (5.28) we have 

E E ^ ^^t!;'^^""' A-O^- ^ J = 2;7^ (5.29) 

and 

(5.30) 
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Hence by using (5.30) and lemma 5.2 it is easy to show the system of equations (5.29) is the 
part of (5.19). Thus here we get only system of equations (5.30) which is (5.12). 

Now consider the last equation, i.e. the equation (4.3). Difference of the left and the right 
sides of equation (4.3) we denote by K, i.e. 

^ " ^ ^^^^ 2(2m-2)! " 2(2m-3)! ^ ^ ~ 2(2m-l)! 

and therefore we keep in mind that K — Q. 

Applying binomial formula for first expression of K and taking into account (4.4)-(4.6), after 
some simplifications 

2(2m-l)! 2(2m-2)! 2(2m-3)!V2 / 2jl{2m - 1 - j)\ 



^ (-1)^- _ (-lyB 



2(2m - - 1)! 2j!(2m - 3 - j)\ 



^2(2m-2-i)!(i-l)!^^(i-l)!(2m-l-i)!:^ ^^^^ 



;^(2m-l-j)!(j-l)! 



(n, - 1 V+1 9m 



. (5.31) 



N 

Now consider the sum ^ C['y]{h'-fY'^'~^''^ in (5.31). For this sum using formulas (5.1), (2.4)- 

7=0 

(2.6) after some simplifications we obtain 

^ ^ ' U ^!(2m-j-i)! 

_ ^T^- (2m-l-i)!/.^^^ ^ ^ d,q--^+p,q,{-ir\ 

^ p-(2m-l-i-p)!^^ {l-q,r^ ^'-^^^ 
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Substituting (5.32) into (5.31) we get polynomial of degree 2m — 1 with respect to h. It is easy 
to see that constant term and coefficients in front of h and h'^ are zero. Then for K we obtain 

^(i-l)!(2m-l-i)! 1^ ^ i\{2m-j-i)\ 

^ p!(2m-l-j-p)!^^ {l-quf^' ) 

, V ^ / f _1 V V ^^^k+PkQk (-1) . ^ 2m-l-i I 

^^(i-l)!(2m-l-i)! ^ ^ ^ (l-g.)^+^ 



..1 ..0 (1-^^)^^^ y ^(2m-i)!(i-l)!- 

Hence if take into account lemma 5.2 and (5.4*), then the expression in the second parenthesis 
are simplified and K has the form 



rO'-l)' ^!(2m-j-z)! ^(2m-j)!(j-l)! 



m— 1 



1. p!(2m-l-i-p)!^- ^'••''^ 

Using lemma 5.3 from (5.33) we get 



j\ ^ i\ {2m — 1 — j — i)\ ^ i\ ^ i\{2m — \ — i — i)\ 

j=3 1=0 ^ J I j=m+l j=0 ^ J I 

m-l 7 2m-j n rn+l , „ m-2 . 

(9.m. - ill - m ^ 2^ - m 2^ 



^ (2m - j)! (j - 1)! (j - 1)! j^^t\{2m-l-j- ^)! 

.ir+2 0'-l)' Z. ^l{2m-l-J-^)\ 
Hence using (5.19) after simplifications we have 



{m — 1)\ \ m J i\ (m — 1 — i)\ 
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h^Z,., (2m - 1 - m-iy h'"^-'Z2m-2 



'I'-^j-i 

^ a -l)\(2m-l- i)\ ^ 



,=n.+2(^'-l)'(2"^-l-^')! U ^!(2m-l-J-^)! (2m -2)! 



[(m-l)!f U ^"^"V .._4l,i!(2m-l-j)!^^ ^ ^ (2m-: 



2m-2 

[(m - 1)!]' V m "^""V ' j! (2m - 1 - j)! ^' ^' (2m - 2)! 

Here the middle sum is equal to zero because when j is even (— 1)^ — 1 = 0, and when j is odd 
Bernoulli numbers = 0. Therefore finally for K we get 

(2m -2)! +[(m-l)!]H"^ ^m-ij-^- 
It means hence taking into account (5.4*) we get following equation 

m-l m-l , TV+j I ^ ^ / R 

y y duq, +Pkq,{-1) ^.Q^-i ^ 

m-l 2m-2 , /v+j , ^ ^ / i 
ifc=l j=0 ^ ^'^^ 

Using (5.30) when a — 1 and taking into account lemma 5.2 it is easy to show , that the 
equation (5.34) coincide by equation of the system (5.19) when a—1. Thus, we have obtained 
the last equation (5.35) for and p^, which is the same as (5.13). Theorem is proved. 

5.2 Calculation of norm of the error functional 

For square of norm of the error functional (1.2) of optimal quadrature formulas of the form 
(1.1) take placed following theorem. 

Theorem 5.2. For square of norm of the error functional (1.2) of optimal quadrature 
formula of the form (1.1) following is valid 



4™)*(o,i) =(-1) 



2 



\m+l 



B2n.h'- (l-gf)(4g^ + (-l)'P.gfe) .,„2^ 

(2m)! ^ (2m)! (1 - g.)^+^ 



where dk,Pk are determined from system (5.10)-(5.13), B2m o-fe Bernoulli numbers, qk are roots 
of Euler polynomial of degree 2m — 2, \qk\ < 1. 
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Proof. Computing defined integrals in the expression we get 



= (-1) 



m+l 



A-B 



A-B 



N 



N 



N 



N 



2(2m-2)! 

|/i/3-/i7p"*-i 



+ 5^ C[P] Fm + 5] C[/3] F{hp) - ^ C[7] ' 2(2^ _ IV 

/3=0 /3=0 I 7=0 ^ ^' 



+ 



A -B{hp- 1) 

2(2m-2)! 



2m-2 



(2m + 1)! 



where F{hp) is determine by formula (5.16). As is obvious from here according to (4.1) the 



m—l 



expression into curly brackets is equal to the polynomial P^_i(/i/9) = ^ \a{h(5)°'. Then \ \^\\ 

a=0 

have the form 



(-1) 



m+l 



A-B 



N 



A-B ^ ^ ^^^j A jhpf^-' -Bjhp- If^-^ ^ 



(2m -3)! (2m -1)! ' ^ 2(2m-2)! 



N N ^ 



/3=0 /3=0 

Hence using (4.2) and (4.3) we get 

A-B A-B 



m+l 



(2m -3)! (2m -1)! 



+ A- 



( 1 

"''•I 2(2m-l)l ~^"^" + 



2(2m- 1)! 



2(2m-3)! 



+ Ai- 



+ 



B 



2(2m-3)! 



N N ^ 

+ J] C[(5\ F{hp) + J] p^^.m - 



/3=0 /3=0 

Prom (5.35) after simplifications using (5.17), (5.12) we have 



(5.35) 



(-1) 



■m+l 



B-A 



2(2m-l)! (2m -2)! I 2(2m 



1 1 ^ 

^ + ^E^w(/^7)^"^-^ + 

' 7=0 



+ 



(2m-2)S 



m— 1 



+ i..E 



^ (2m - 1 - j)!(j - 1)! V2(2m - j) 2m - j 

m-l2m-l-j ^ N+i(_T\i -r N 

_h2m-jJ2 J2 ' A^02m-l-i^£^^[^](_;^^)2m-l-j_ 

(?fe 1)* 



fc=l i=0 



7=0 
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(2m - 1 - J)^^ 1 + Yl ^[^] ^('^^) + E ^t^] 



(2m + 1)! 



/3=0 /3=0 

Hence, taking into account (5.16), (5.17), (5.18) and using (4.4)-(4.6), after some calculations 
we obtain 

m— 3 7-1/ -, \ 1 m— 1 



-1) 



m+1 



1 



B i-iy 



+ E77 



(2m -1)! (2m +1)! ^ (2m - 2 - j)! (j + 1)! ^ (j + 1)! (2m - j)! 



m— 1 



-E 

j=2 



m-1 2m- 1 -J 



(2m - 1 - j)!(j + 1)! \ 2m-j 



k=l i=0 



m—1 



E 



AT 



(g;t - 1)^+^ 



^iQ2m-l-j ^ 



-(2m-l-j)!(j + l)!^ 



J]C[7](/i7)'™-^-^+E^W 



/3=0 



Ml 

(2m)! 



(5.36) 



When q; > m — 1 using lemma 5.1 and formulas (2.4)-(2.6) we get 



N a-1 . 



m— 1 a 



fe=l 1=0 



(1 - qkY^' 



^i!(a-j)!A.2.- 



(1 - qkY"-' 



(5.37) 



Using (5.37) and taking into account (2.1), (2.3), after simplications, from (5.36) we have 



+ 



-l)'"+i 



(2m 



■2m-3 

E 

.3=m 

EE 



i?2m-j^^™ , B2mh 



2m 



+ 



(2m- j)!(j + 1)! (2m)! 



' fc=l i=0 

2m- 1 / 



(1 - qkY^^ 

a-2 



.^jQ2m^ 



+ E 



1)" g 



^ (2m-Q;)!^ 7!(a + l-7')! 

a=m ^ ' 1=1 ^ •'^ 



2m 



-E 



m-1 i 

(2m-a)!^j!(a-j)! 



■1)" 



49r^+Pfcg;c(- l) 
(1 - 5;.)^+^ 



i+l 
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Hence last two sums regrouping in powers of h, using designation (5.4*) and keepin in mind 
Za-i = a = 3, m, we have 



m+l 



■2m-3 



L,. . ... + + 



(2m-i)!(j + 1)! (2m)! 



EE 



2m- 1 



a! ^ (2m -?■)!(?■- a + 1)! ^ a\ ^ (2m - 7)! (7 - a + 1)! 



. X- 2m+l „ 2m 

('9m - I'll ('-i - -U I'll Z^ ('n/ - I'll 2-^ 



, a! ^ (2m- 7)! (7 - a + 1)! ^ fa-l)! ^ (2m - 7)! (7 - a + 1)! 

a=3 j=m ^ J J \J 'J a=m+l ^ 3=a-l ^ J J \J 'J 



(5.38) 



Since in (5.38) 

2m 



E 



(-1)^ 



a-l 



-^_^{2m- - a+l)\ (2m-a-l)! 



(1-1) 



2m-a+l 



0, m + 1 < q; < 2m 



and 



2m- 1 

E 

j=a 



[-ly 



(2m - j)\ {j -a + l)\ (2m - a + 1)! 



((-ir^ + 1), 



then, using these equalities, from (5.38) we obtain 



:-i) 



m+l 



2J /^2m ^2m+l "^-^ 2m 



+ 



(2m)! (2m)! ^ 



EE 



i=0 



(l-gf)(49^ + (^l)-pfegfe) 



2m-l 



_ B^ ii-ir-' + 1) 

Z^ a\ (2m- q;+ 1)! 

jr=m+l ^ ' 

Hence taking into account that when a is even (—1)"^^ + 1 = and when a is odd B^ = 0, 
(since a ^ 1), we get the statement of theorem 5.2. Theorem 5.2 is proved. 

Corollciry 5.1. In the space L2 (0,1) among quadrature formulas of the form (1.1) with 
the error functional (1-2) there exists unique optimal formula which coefficients are determined 
by following formulas 

C[I3] - ' ^ 



h, (3=1,N-1, 
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A=—, B^- — . 
12 12 

Furthermore for square of norm of the error functional following is valid 

2 h"^ 



720' 



CoroUciry 5.2. In the space L2^^(0, 1) among quadrature formulas of the form (1-1) with 
the error functional (1-2) there exists unique optimal formula which coefficients are determined 
by following formulas 



h 

2 ' 


f3 = 


0,N, 




= 


1,N-1, 


12' 


B = 


12' 



A 

Furthermore for square of norm of the error functional following is valid 



£|Lr(0,l) 



30240 

Proofs of Corollaries 5.1 and 5.2 we get immediately from theorems 5.1 and 5.2 when m = 2 
and m = 3 respectively. 
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